A multi-wavelength approach to the properties of Extremely Red 
Galaxy populations: I - Contribution to the Star Formation Rate 

density and AGN content. 
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ABSTRACT 

We present a multi- wavelength analysis of the properties of Extremely Red Galaxy (ERG) popu- 
lations, selected in the GOODS-South/Chandra Deep Field South field. By using all the photo- 
metric and spectroscopic information available on large deep samples of Extremely Red Objects 
(EROs, 645 sources), IRAC EROs (IEROs, 294 sources), and Distant Red Galaxies (DRGs, 350 
sources), we derive redshift distributions, identify AGN powered and Star- formation powered 
galaxies, and, using the radio observations of this field, estimate robust (AGN- and dust-unbiased) 
Star Formation Rate Densities (p*) for these populations. We also investigate the properties of 
"pure" (galaxies that conform to only one of the three ERG criteria considered) and "combined" 
(galaxies that verify all three criteria) sub-populations. Overall, a large number of AGN are iden- 
tified (up to ~ 30%, based on X-rays, and mid-infrared criteria), the majority of which are type-2 
(obscured) objects. Among ERGs with no evidence for AGN activity, we identify sub-populations 
covering a wide range of average star-formation rates, from below 10 Mq yr _1 to as high as 200 
Mq yr _1 . Applying a redshift separation (1 < z < 2 and 2 < z < 3) we find significant evolution 
(an increase of a factor of 2 or higher) of for EROs and DRGs, while none is observed for 
IEROs. The former populations can contribute more than 20% to the global p» at 2 < z < 3. The 
emission from AGN activity is typically not strong in the ERG population, with AGN increasing 
the average radio luminosity of ERG sub-populations by, nominally, less than 20%. AGN are 
common, however, and, if no discrimination is attempted, this could significantly increase the p* 
estimate (by over 100% in some cases). Thus, and while the contribution of star forming processes 
to the radio luminosity in galaxies with AGN remains uncertain, a comprehensive identification 
of AGN in these populations is necessary to obtain meaningful results. 

Subject headings: galaxies: high-redshift; galaxies: active; galaxies: starburst 
1. Introduction 
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The study of high redshift galaxy populations is 
fundamental in understanding the formation and 
evolution of galaxies. One of the major difficul- 
ties in such studies is assembling an unbiased and 
representative sample of galaxies at high-redshift. 
With the advent of highly sensitive infrared detec- 
tors, both ground- and space-based, it has become 
clear that optical observations reveal only a small 
and very biased fraction of the galaxy population. 
This started with the Infrared Astronomical Satel- 
lite (IRAS) in the 1980s, which found a population 
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of highly luminous galaxies, radiating most of their 
immense luminosity at infrared w avelengths (1- 
lOOOjLtm. lSanders fc Mirabel Ill996l ). At about the 
same time, ground-based telescopes using increas- 
ingly more efficient near-infrared (NIR, 1-2.2/im) 
detecto rs, allowed the detection of very red galax - 
ies (e.g. lElston et al.lll988tlThompsori et~aDll999h . 
More recently, observations at (sub-) millimeter 
wavelengths revealed a population of distant mas- 
sive galaxies, unseen at optical wavelengths, go- 
ing through intense (and likely short) episodes 
of star- formation (SF) activity, possibly leading 
to the formation of present-day ellip ticals (e.g. 
Barger et al1ll998t iHughes et al.lll998l ). It is now 
clear that understanding the z > 1 universe re- 
quires deep multiwavelength observations to build 
up a complete picture of the processes occurring 
early in the evolution of galaxies. 

In an attempt to constrain hierarchical models 
of galaxy formation, the last few years have seen 
optical-to- infrared or infrarcd-to-infrared color cri- 
teria being used to find high-redshift galaxies 

hosting evolved stell ar population s. Extremely 

Red Objects (EROs. iRoche et al-lEx)! . IRAC- 
selected EROs (IERO s, also known as IR Ex- 
tremely Red Objects, lYan et al.l 12004) and Dis- 
tant Red Galaxies (DRGs. iFranx et al.l 120031 ) 



were thought to identify old passively evolving 
galaxies at increasing redshifts (from z > 1 for 
the EROs/IEROs to z > 2 for the DRGs), for 
which a prominent 4000 A break would fall be- 
tween the observed bands. These techniques, 
however, are also sensitive to active (star-forming, 
AGN or both) high-redshift dust-obscured galax- 



butions dSm< 


fil et al. 2002; Alexander et al. 2002; 


Afonso et al. 


2003; Papovichet al. 20061). These 



active members of the Extremely Red Galaxy pop- 
ulation (ERGs, as we will collectively call EROs, 
IEROs, and DRGs) are also important targets for 
further study given that they constitute a dusty 
population of g alaxies easily miss ed at optical 
wavelengths (e.g. Afonso et al.l l2003). A challenge 



in studying the nature of the red galaxy popula- 
tion is the difficulty in disentangling the effects 
due to redshift, dust-obscuration, and old stellar 
populations. 

Identifying the so-called "Passively Evolving" 
and "Dusty" ERGs is a fundamental and partic- 
ularly difficult task, where optical spectroscopic 



observations are of limited use. The identifica- 
tion and study of Active Galactic Nuclei (AGN) 
or star-formation (SF) activity in these galaxies, 
for example, requires multi-wavelength data from 
X-ray to radio wavelengths. Radio observations 
are of particular interest here, given the possibil- 
ity to reveal the activity in these obscured systems 
and, for star-forming dominated galaxies, allowing 
for a dust-free estimate of their SF rates (SFRs). 

In this paper we present a comparative study of 
the ERG population. Using the broad and deep 
wavelength coverage in the Great Observatories 
Origin Deep Survey (GOODS)-South/Chandra 
Deep Field South field, we select samples of EROs, 
IEROs, and DRGs and estimate their statisti- 
cal properties. With the extensive photometric 
data available we explore the redshift distribu- 
tion, SFRs, and AGN activity in these galaxies. 
Using radio stacking we estimate dust-free SFRs 
and the contribution of red galaxy populations 
with no detected AGN activity to the global SFR 
density (/}*). 

The structure of the paper is as follows. Sample 
selection is described in Section [2] Section [3] ad- 
dresses the AGN identification technique. In Sec- 
tion [4] the ERG sample is characterized, leading 
to the estimate of the dust- and AGN-unbiased 
contribution to p*. Finally, the conclusions are 
presented in Section [5] 

Throughout this paper we use the AB mag- 
nitude system^, and a ACDM cosmology is as- 
sumed with Ho = 70 km s _1 Mpc -1 , CIm = 0.3, 
0\ = 0.7. 

2. Sample Selection 

GOODS was designed to assemble deep multi- 
wavelength data in two widely separated fields: 
the Hubble Deep Field North (HDF-N/GOODS- 
N) and the Chandra Deep Field South (CDF- 
S/GOODS-S). Specifically the southern field - 
GOODS-S — includes X-ray observations with 
Chandra and XMM-Newton; optical (BVIz) high 
resolution imaging with the ACS on board the 



^^When necessary the following relations are used: 
(K, H, J, I ) AB = (K, H J, I ) Vega + (1.841, 1.373, 0.904, 
0.403) from lRoche et al.l j2003l ); 

IRAC: ([3.6], [4.5], [5.8], [8.0]) AB = ([3.6], [4.5], 
[5.8], [8.0]) Voga + (2.79, 3.26, 3.73, 4.40) from 
http://spider.ipac. caltech. edu/staff / gillian/ cal.html 
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Hubble Space Telescope HST; near-infrared (NIR) 
and mid-infrared (MIR) coverage with the Very 
Large Telescope (VLT) and the Spitzer Space 
Telescope, respectively; and radio imaging with 
the Australia Telescope Compact Array (ATCA), 
Very Large Array (VLA), and Giant Metrewave 
Radio Telescope (GMRT). These data are among 
the deepest ever obtained. Large programs aim- 
ing at comprehensive spectroscopic coverage of 
this field are also being performed. The quality 
and depth of such data make these fields ideal to 
perform comprehensive studies of distant galaxies 
and, in particular, of the ERG population. 

We use the FIREWORKS X s -band selected 



2.1. Red Galaxy Samples 

We consider three categories of ERGs: 



catalog from GOODS-South (jWuvts et al.ll2008ft . 
This provides reliable photometry from UV to 
IR wavelengths (0.2-24/xm) for each source de- 
tected in the K s ISAAC/VLT maps (ISAAC 
GOODS/ADP vl.5 Release), thus covering an 
area broadly overlapping the HST ACS observa- 
tions. The widely different resolutions between 
optical and infrared bands are properly handled 
to allow consistent color measurements. This is 
performed by adjusting the optical HST and NIR 
VLT images to a common resolution, and per- 
forming photometry on optical, IRAC, and MIPS 
images, using the prior knowledge about position 
and extent of sources from the K s -ha,\\A images. 

Redshift estimates are also provided. We only 
use spectroscopic observations to get redshifts 
with which to derive the intrinsic luminosities of 
ERGs. Only good spectroscopic redshift deter- 
m inations (quality fl ag equal or greater than 0.5 
in lWuvts et al.l[2008h were considered, comprising 
21% of the ERGs (Table HJ. For the remaining 
sources, photometric redshift estimates from the 
FIREWORKS catalog were considered. The red- 
shift distributions will be discussed in Section HtI 



The FIREWORKS catalog contains 6308 K s 
selected sources. To allow for robust selection of 
our ERG populations we impose a magnitude com- 
pleteness limit of K s ^ T o T = 23.8 AB (S/N> 3a) 
and, following the p rescription for robu st photo- 
metric samples from Wuvts et al. (2008), adopt a 
pixel weight limit of K s w> O.sH. This results in 
a final ^-selected catalog of 4740 sources with 
robust K s photometry. 



• EROs: j 775 -K s > 2.48 (|Roche et al-llioO 



• IERO s: z 850 - [3.6] /jm > 3.25 (|Yan et al 



2004); 



DRGs: J — K s > 1.30 (jFranx et alj |2003). 



2 See Section 3.4 of IWuvts et al.l J200STI for a description of 
the concept of pixel weight. 



Whenever a source is not detected in one of 
the bands, a limit to its magnitude is assumed 
(we adopt the 3a flux level based on the local rms 
provided in the catalog). In the case of unreliable 
photometry (e.g., K s w< 0.3), the corresponding 
source is not considered further. Thus, only ERGs 
with robust photometry in both bands (or with 
a non-detection in the shorter wavelength band) 
used for their identification are considered. 

The resulting ERG sample, with robust pho- 
tometry, contains 731 objects: 645 EROs, 294 
IEROs, and 350 DRGs, down to the adopted mag- 
nitude limit of K St TOT — 23.8 AB. These classi- 
fications overlap, with individual objects poten- 
tially included in more than one classification. 
This is illustrated in Figure Q] 

It should be noted that FIREWORKS is a K s 
selected catalog. As such, EROs and DRGs are se- 
lected according to the traditional definition, but 
IEROs selected from the FIREWORKS catalog 
are in effect if s -detected IEROs. This sample will 
only be representative of the true IERO popula- 
tion in the absence of a significant number of very 
red K s — [3.6] IEROs, which are undetected in the 
K s image. One should also note that the zsso de- 
tection limit (the current K s selected sample in- 
cludes sources with up to z^q ~27mag) imposes a 
[3.6]-band magnitude limit of ^23.75 mag for the 
IERO sample. 

Figure [5] shows the color-magnitude distribu- 
tion for sources in the FIREWORKS catalog and 
for the ir s -detected IERO sample, displaying our 
adopted i<T s -band magnitude limit (diagonal line) 
and the practical [3.6]-band magnitude limit (ver- 
tical line). The sampled region at if^TOT — [3-6] < 
-[3.6] +23.8 and [3.6] < 23.26 (below the diagonal 
line and to the left of the vertical one) does not 
indicate a major incompleteness for the critical re- 
gion (above the diagonal line and to the left of 
the vertical one). For example, allowing all FIRE- 
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WORKS sources to be considered (up to a K s - 
band magnitude of 24.3 mag), would only increase 
the IERO sample by 6% (17 new IEROs). Conse- 
quently, we consider our sample of (if s -detected) 
IEROs representative of the true IERO popula- 
tion and find it unnecessary to assemble a sepa- 
rate sample of IEROs from a 3.6 /im selected cata- 
log, thus maintaining the photometric homogene- 
ity within the ERG sample. 

2.2. Sub-classes of ERGs 

We refer to those sources that appear in only 
one of these classes (ERO, IERO, or DRG) as 
"pure" populations, while those that are simulta- 
neously included in three ERG categories are re- 
ferred to as the "common" population. In this 
work, the latter will be referred to as "common" 
ERGs, or cERGs. When addressing both the 
"pure" and "common" populations we will restrict 
ourselves to those sources which have sufficient 
information for a definitive classification in each 
of the three red galaxy criteria (either good pho- 
tometry or robust upper limits in all bands used 
for classification). With such requirements, we 
find 702 ERGs: 626 EROs, 289 IEROs, and 338 
DRGs. Essentially all IEROs (290 out of 29$) 
are also EROs and more than two-thirds (198 out 
of 289) are also classified as DRGs; there are 267 
sources that comply with both the ERO and DRG 
criteria, and 197 ERGs that are simultaneously 
classified as ERO, IERO, and DRG (the cERGs). 
Identified as "pure" sources are 272 pure EROs 
(pEROs), 3 pure IEROs (pIEROs), and 70 pure 
DRGs (pDRGs). 

Figure[T]shows the overlap between the different 
sub-populations. The initial columns of Table [T] 
summarize the numbers referred to above. 

3. Mult i- wavelength AGN identification 
and classification 

One of the major problems for the characteriza- 
tion of ERGs, or for any distant galaxy population, 
is to identify the presence of AGN activity. The 
many techniques that exist target different AGN 
types and redshift ranges, and no single technique 
can guarantee a high discriminatory success rate. 



3 Number of IEROs with good photometry in the 1775, 2850 > 
K s and 3.6/^m bands. 



X-rays, radio or MIR, originating from different 
regions in the vicinity of the AGN, and differently 
affected by dust obscuration, provide independent 
ways to reveal such activity. Hence, we use the 
multi-wavelength data available in this field to 
carry out a thorough identification and classifica- 
tion of AGN activity in the ERG population. 

3.1. Optical Spectroscopy 

Optical line ratios will only reveal AGN activ- 
ity if most of the galaxy's line emission comes 
from the environment near the AGN and if the 
dust obscuration is not significant. In the case 
of obscured AGN activity, the emission from any 
disk SF may dominate the optical line emission. 
Also, since ERGs are intrinsically UV/optically 
faint, spectroscopy will be of limited use to re- 
veal their nature. Overall, there are only 11 spec- 
troscopic AGN identifications (narrow line AGN 
or QSO type-2 classifications), galaxies which are 
also identified as AGN by the criteria described in 
the following sections. 

Spectroscopy also allows for the rejection of 
galactic stars selected as EROs. In this sample, 
4 were found and discarded from further study. 

3.2. X-Rays 

X-ray emission is arguably the most effective 
discriminator of AGN activity in a galaxy. Due 
to the sensitivity levels currently reached with 
the deepest observations (the 2 Ms C DF fields: 
Alexander etaD 120031 iLuo et al.l Eooih . the most 
powerful AGN (£o.5-8fceV > 10 43 ergs _1 ) can be 
detected beyond the highest redshift currently ob- 
served, z > 7. On the other hand, both low lu- 
minosity AGN and vigorous starforming galaxies 
(£o.5-8/c e v ~ 10 41_42 erg s^ 1 ) can only be de- 
tected out to z ~ 1 — 2. If enough signal is de- 
tected, detailed spectral analysis can be used to 
distinguish between AGN and SF activity as the 
origin of the X-ray emission. 

In this work, the ERG sample was cross- 
matched with the rece nt catalogs from the 2 Ms 
Chandra observations ( Luo et al. 20081) . For the 
region considered here - GOODS-S ISAAC - the 
X-ray observations reach aim-point sensitivity lim- 
its of« 1.9x 10~ 17 and « 1.3 x 10" 16 ergc m - 2 s" 1 
for the soft (0.5-2.0 keV) and hard (2-8 keV) 
bands, respectively. 
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X-ray detections were searched for within 2" of 
each ERG position. Counterparts were found for 
84 of the 645 EROs (-13%), 47 of the 294 IEROs 
(-16%), and 51 of the 350 DRGs (-15%) (see 
Tabled]). These dete ction fractions ar e consi stent 
with those found by lAlexander et al.l ([2002), for 
EROs, and lPapovich et al.l ([20061 ). for DRGs. 

We adopt the X-r ay classification criteria from 



Szokolv et al.l (|2004h . which considers both the 
X-ray Luminosity (Lx), estimated from the 0.5- 
8 keV flux, and hardness ratio (HR), calculated us- 
ing the count rates in the hard band (HB, 2-8 keV) 
and in the soft band (SB, 0.5-2 keV) bands: HR 
= (HB-SB)/(HB+SB). This is listed as follows: 

Galaxy : L x < lO^ergs" 1 & HR < -0.2 
AGN - 2 : 10 41 < L x < lO^ergs" 1 & HR > 
AGN - 1 : 10 42 < L x < lO^ergs" 1 & HR < 
QSO - 2 : L x > lO^ergs^ 1 & HR > -0.2 
QSO - 1 : L x > lO^ergs" 1 & HR < -0.2 



0.2 
0.2 



The rest-frame X-ray luminosity is calculated 



as: 



L 



x 



47r4/ x (l + z) r - 2 ergs- 



where dj, is the luminosity distance in (cm), fx 
is the X-ray flux (erg s^ 1 cm~ 2 ) in the 0.5-8 
keV band and the photon i ndex is assumed to be 



r = 1.8 dTozzi et al.l [20061. The luminosity dis- 



tance is calculated using either the spectroscopic 
redshift or, if not available, the photometric red- 
shift estimate. 

In total, these criteria enable the identification 
of 86 sources hosting an AGN with only 6 X-ray 
sources powerful enough to be classified as QSOs. 
The majority of the AGN are classified as type- 
2 sources: 41 X-ray detections have HR > —0.2 
(i.e., obscured) while only 19 show lower values 
(with the remaining 26 having uncertain HR de- 
terminations, with no discrimination possible), in- 
dicating a possible 2:1 obscured to unobscured 
ratio. However, no constraints on the ratio of 
type-2 to type-1 AGN can be established, since 
HR becomes degenerate as an obscuration mea- 
sure above z — 2, except for the most heavily ob- 
scured AGN(A^h ^> 1 23 cm~ 2 ; see, for example, 



Alexander et al.ll2005l ). Figure [3] shows the HR- 
redshift evolution for simple obscured and unob- 
scured AGN X-ray emission models, together with 
the X-ray detected AGN ERGs. Although Type-2 
sources seem to dominate (Table [T]), an X-ray es- 
timate of the Type-2 to Type-1 AGN ratio among 



such high-redshift populations would require ob- 
servations extending to lower X-ray energies (<0.5 
keV) below those reliably achieved by Chandra. 

3.3. Mid-Infrared 

Over the last few years with the sensitiv- 
ity of IRAC and MIPS onboard Spitzer, several 
MIR criteria have been developed for the iden- 
tification of AGN at the center of galaxies. A 
power-law (PL) MIR spectral energy distribution 
(SED), for example, is charact eristic of AGN emis- 



sion (e.g. iDonlev et al.l 120071) . Somewhat more 



generic color-color diagrams have also been inves- 
tigated, and AGN loci in such plots defined (e.g 
| lvison et all 12004 lLacv et al.l 12004 IStern et al 



2005; lHatziminaoglou et al.l 120051) . This wave- 
length range is of particular interest for the ERG 
population, given their red SEDs. Here, we have 
applied MIR diagnostics to our ERG sample, as 
described below. 

Observational data at X-ray and IR wave- 
lengths provide complementary views of AGN ac- 
tivity. The most obscured AGN may be missed 
by even the deepest X-ray surveys but can still be 
identified by their hot-dust emission at IR wave- 
lengths. On the other hand, depending on the 
amount of dust and its distribution, and on the 
AGN strength, the MIR emission from X-ray clas- 
sified AGN may not be dominated by the hot dust 
in the vicinity of the AGN itself. A detailed com- 
parison of the merits of AGN selecti on by the X- 



rays a nd the MIR was performed bv lEckart et al 



(l2010f ). showing that only a multiwavelength com- 
bination of AGN criteria can help to overcome 
biases present in single-band selection. However, 
even the combination of MIR and X-rays will not 
result in complete AGN samples, as the identifi- 
cation of low power AGN will ultimately depend 
on the depth of the surveys. By performing this 
study in GOODS-S, with some of the deepest data 
both at X-ray and MIR wavelengths, we maximize 
the identification rate of AGN. 

IRAC counterparts were found for practically 
all (98%) ERG sources. The vast majority (92%) 
are detected simultaneously in all IRAC bands: 
612 of the 645 EROs, 290 of the 294 IEROs, and 
322 of the 350 DRGs. The MIPS 24^m detec- 
tion rate is understandably lower (55%), given 
the lower relative sensitivity: 359/191/222 of the 
645/294/350 EROs/IEROs/DRGs are detected 
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(Table EE}. 

3.3.1. Classification: MIR spectral index 

The IRAC fluxes for each ERG, covering the 
3.6 - 8.0 /jm range, were fitted with a PL (f u oa ,a ), 
classifying as AGN sources those where a < —0.5 
rtDonlev et al To increase its reliability, ap- 

plication of this technique was limited to sources 
simultaneously detected in all of the four IRAC 
bands, and the a value was only accepted if the 
X 2 probability fit was P x 2 > 0.1. 

This criterion reveals the existence of AGN in 6 
(1%) EROs, 1 of which is also identified as an AGN 
in X-rays; in 4 (1%) of the IEROs, 1 of which have 
an X-ray AGN classification; and in 12 (3%) of the 
DRGs, none having an X-ray AGN classification 
(see Table EE}. 

3.3.2. Classification: MIR colors 

In recent years, several AGN color-selection cri- 
teria have be en developed employing MIR IRAC 
observation s ( Ivison et al. 1 12004 lLacv et alj|2004 
IStern et all 120051: lHatziminaoglou et all 12005^ . 
They are not completely independent from the 
PL criteria mentioned above, but are likely more 
sensitive to other kinds of AGN (e.g., QSOs or 



broad-line AGN, lLacv et al 



Stern et al.l 

2005) . Her e we fo llow the criterion proposed by 



2004t [ 



Stern et al.l (|2005D . which identifies AGN from 
galaxy populations if they satisfy the following 
relations: 



([5.8] - [8.0]) > -0.07 
([3.6] - [4.5]) > 0.2([5.8] 
([3.6] - [4.5]) > 2.5Q5.8] 



[8.0]) - 0.156 
[8.0]) - 2.295 



Figure 2] shows t he distributi o n of E RGs on the 

" (|2005l) . It identi- 



color-color plot of IStern et al 



fies as AGN 135 (21%) EROs, 27 of which are 
also classified as AGN from the X-rays; 83 (28%) 
IEROs, 21 of which also have an X-ray AGN clas- 
sification; and 117 (33%) DRGs, 25 of which also 
appear as X-ray AGN. The relatively high number 
of potential AGN identified, over th at revealed by 



the X -rays, is known and expected (|Donlev et al 
20071 and references therein). 



3.3.3. MIR degeneracy at z > 2.5 

One problem in using MIR photometry to iden- 
tify AGN (with both PL and color-color criteria) 



arises at z > 2.5, as both star- forming galaxies and 
AGN start to merge into the same MIR color-color 
space. The main reason for this is the increasing 
relative strength of stellar emission in the MIR, 
as compared to that of an AGN, as redshift in- 
creases. At higher redshifts, a prominent 1.6 (im 
stellar bump may pass through the IRAC bands, 
allowing for the detection of a steep spectral in- 
dex not from AGN emission, but from the stellar 
emission alone. 

At z>2.5, the MIR spectral index method (see 
SectionEXI} identifies AGN in 3 EROs, 3 IEROs, 
and 10 DRGs, where none is identified through X- 
ray emissio n; on the other ha nd, the color-color 
method of IStern et al. (|2005l) classifies as AGN 
51 EROs, 42 IEROs, and 76 DRGs, with 14, 11, 
and 16 (respectively) X-ray confirmed at these red- 
shifts. 

In the present work, this is not a serious prob- 
lem, as most of the ERG sample (77%) lies at 
z < 2.5 (see Section l4~Tj) . Nevertheless, it should 
be noted that at higher redshifts, this could result 
in a likely overestimate of the presence of AGN. 
One can attempt to correct for this effect, by us- 
ing the MIPS 24/im observations: at z ~ 2.5 — 5, 
the 1.6 /mi bump will be shifted to the 6-10 /im 
range. Therefore, in the absence of significant 
AGN emission, one expects a blue [8.0]- [24] color. 
We have investigated the [8.0]- [24] colors of z > 2.5 
ERGs, removing a MIR AGN classification when- 
ever [8.0]- [24] < 1. In this case, the observed MIR 
emission does not require an AGN component, but 
just a more or less evolved stellar population. 

In Figure Owe present the MIR [5.8]-[8.0] ver- 
sus [8.0]-[24] color-color plot for z > 2.5 ERGs 
in the current sample (highlighting those classi- 
fied as AGN by the two MIR criteria referred to 
in Section 13.3.11 and Section I3.3.2[) . The tracks 
represent the expected colors of template SEDs 
when redshifted between z ~ 2.5 and z = 4. 
All templates were taken from the SWIRE Tem- 
plate Library (jPolletta et al.l 120071 ). with the ex- 



ception of the track in the upper right, corre 



spond ing to the extreme ERO of lAfonso et al 



(|200ll ), the MIR emission of which is dominated 
by an obscured AGN. The vertical li ne indicates 



the [5 .8]-[8.0] color constraint of the IStern et al 



(|2005f ) criterion, while the horizontal line shows 
our adopted color cut separating AGN and star- 
forming processes at these redshifts. The AGN 
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template that crosses over this [8.0]- [24] thresh- 
old at the highest redshifts is IRAS 22491-1808, a 
mixture of AGN and stellar MIR emission, where 
the AGN component is progressively less sampled 
by the MIR bands as redshift increases. Con- 
cerning the current sample, the few high redshift 
(z > 2.5) sources classified as AGN by the MIR 
criteria that appear below the [8.0]-[24] thresh- 
old, 10/7/10 EROs/IEROs/DRGs, do not require 
an AGN SED to explain their MIR emission and, 
consequently, their MIR AGN classification is re- 
moved. We note that the MIPS (24/zm) - IRAC 
colors have not yet been considered in detail as a 
way of selecting AG N/SF sources in deep MIR sur- 
veys (however, see llvison et al. 2004; lLacv et al.l 
2004h . Prompted by Figure [5] we are currently 
performing a detailed analysis of such colors as 
potential diagnostics for AGN behavior at high-z, 
to be presented in a subsequent paper. 

We note the presence of two interesting sources 
in Figure [5] The one isolated in the upper right, 
is one of the se ven optically unidentified radio 
sources found in lAfonso et al.l (2006. their source 
#42). Inspection of the K s and 24/im images re- 
veals no signs of blending, strengthening the ac- 
curacy of the 24/zm flux. This source also has 
X-ray emission characteristic of an AGN (Lx = 
10 43.3 ergs -i) The FIREW QRKS assigned photo- 
metric redshift is z = 3.69 (in agreement with our 
own detailed ongoing investigation of this source, 
resulting in an estimate of z — 3.85 using Hyper- 
z). The high-z obscu r ed A GN scenario postu- 
lated 



Afonso et al 



(2006) for this source is 
thus strengthened. The color-track closest to this 
source in Figure [5] is that of the highly obscured 



AGN ERO found by lAfonso et al.l (|200lh 



The other interesting source is the bluest [8.0]- 
[24] 24^m detection, with MIR colors characteris- 
tic of spiral galaxies. It is also X-ray detected but 
has no radio emission. This is a candidate for a 
high-z evolved system with MIR colors typical of 
Spiral c type. The redshift assigned to this source, 
Zphot = 2.53, is at the upper limit of the redshift 
range i n which these type of s ources have ever been 



found ([Stockton et al.l 120081 ). A few more galax 



ies fall in the same region of the color-color plot. 
Again they are candidates for Spiral type evolved 
systems whose study is relevant to constrain hier- 
archical models for the formation of galaxies. 



3.4. Radio 

Radio emission is essentially unaffected by dust 
obscuration, making it a highly desirable diag- 
nostic for SF activity in ERGs. Since both SF 
and AGN activity can produce radio emission, 
it is often difficult or impossible to rely on ra- 
dio properties alone to reveal the power source 
in a galaxy. Indications from radio spectral in- 
dices are of limited use, as both SF and AGN 
emission usually result from synchrotron radiation 
with S u oc ^~ - 8 , and only some AGN show signs 
of flat or even inverted radio spectra. Very high 
resolution VLBI radio imaging has also been used 
with limited success to impose limits on the size of 
the radio emitting region, identifying star-forming 
galaxies where the radio emi ssion is resolved, and 
a possible AGN where not (Muxlow et al.l 2005; 



iMiddelberg et al.ll2008t ISevmour et al.ll2008| ). The 
only straightforward radio AGN criterion is the 
radio luminosity itself, as the highest luminosities 
can only be produced by the most powerful AGN. 

lAfonso etail (|2005h performed a detailed study 
of the sub-mJy radio population, and found star- 
forming galaxies with radio luminosities up to 
ii.4GHz ~ lO^WHz" 1 . We thus take this 
value as the upper limit for SF activity. We note 
that this value correspond s to an SFR of almost 
2000 M Q yr 1 (|Bellll2003l see Section WM- The 
existence of galaxies with higher rates of SF activ- 
ity is unlikely. 

For the current work we have used the 1.4 GHz 
ATCA observations of this field, which reach a 
uniform 14-17 niy rms throughout the GOOD S-S 
field (see lAfonso et al.ll2006l: iNorris et al.ll2006l for 
more details). 

We used as 3"matching radius to search for ra- 
dio emission from the ERGs. In total 17 ERGs 
have reliable radio detections (> 4.5er): 15 EROs 
(2.5% of this population), 10 IEROs (3.3%), and 
11 DRGs (3.1%). Only 3 of these sources have ra- 
dio luminosities in excess of 10 24 5 WHz -1 . They 
are also classified as AGN by the previous X-ray 
and MIR criteria. On the other hand, only 3 radio- 
detected ERGs are not classified as AGN by any of 
the adopted criteria. For these sources, the radio 
luminosity corresponds to SFRs of 93, 646, and 



1584 Mq yr (using the conversion from Bell 



2003. see Section PO 



The small detection rate indicates that power- 
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ful AGN and the most intense starbursts are not 
common in the ERG population, as only sources 
with Li.4GHz > 10 23 WHz _1 will be detected at 
z > 1 with the sensitivity available even in the 
current deepest radio surveys. 

4. Properties of ERGs 
4.1. Redshift Distributions 

As noted above, robust spectroscopic redshifts 
are available for around 21% of the ERG sample. 
Photometric redshift estimates are also available 
from the FIREWORKS catalog, covering almost 
the complete ERG sample. 

The redshift distributions for the ERO, IERO, 
and DRGs are shown in Figure [6] Although the 
range of redshifts sampled in all ERG classes is 
similar, the average value increases from z = 1.76 
for EROs, to z = 2.11 for IEROs, and to z = 2.40 
for DRGs p opulations (in agreement with previous 



work s, e.g. iConselice et al.l 120081 ; iPapovich et al 



2006). This was expected given the source selec- 
tion, designed to identify objects at such redshifts. 

The AGN in the ERG population follow a sim- 
ilar redshift distribution but the AGN fraction in- 
creases rapidly at higher redshifts. This will be 
addressed in the next section. 

Figure [7] displays the redshift distributions once 
again but for pEROs, pDRGs, and cERGs (only 
3 pIEROs exist, at redshifts z p hot = 1.7,3.0 and 
z S pect = 2.6). The redshift distribution of pEROs 
is quite narrow, selecting sources essentially at 
z = 1 — 2 (peaking at z ~ 1.3), while the pDRG 
population is notably small, and at higher red- 
shifts (z = 2 — 4). The "pure" criteria thus ap- 
pear to be good and easy techniques to select 
high-z sources in narrow distinct redshift bins. 
Sources classified as cERGs, appearing as red in 
all three ERG selection criteria, cover a broad red- 
shift range, from z = 1 to z = 4. There are no 
cERGs at z < 1, in this particular sample due to 
the IERO criterion. 

4.2. AGN content of ERGs 

As described in the previous section, several 
multi-wavelength indicators were used to identify 
AGN in the ERG population. The indicators have 
different sensitivities to AGN characteristics, such 
as distance, dust obscuration, or AGN strength. 



Their combination will, thus, allow for a more 
complete census of AGN content in these sources. 

We do not find a numerous population of very 
powerful AGN among the ERGs, as given by the 



X-rays (L x > 10 44 erg s"\ 6 ERGs) and ra- 
dio (Li.4ghz > lO^WHz" 1 , 3 ERGs) luminosi- 
ties. These represent, respectively, only 0.8% and 
0.2% of the ERG sample, comparable to the frac- 
tion observed in the complete -RT-selected FIRE- 
WORKS sample, where 17 (0.4%) QSOs and 5 
(0.1%) radio-powerful sources are found. Overall, 
we select 214 (29%) AGN-dominated systems in 
the ERG sample (29% for EROs, 35% for IEROs, 
and 39% for DRGs). This fraction increases from 
low to high redshift, from 22% at 1 < z < 2 
to 38% at 2 < z < 3. Among the X-ray iden- 
tified AGN, 35% are also classified as such by 
the Stern 'wedge'. Conversely, only 18% of the 
'wedge' identified AGN are X-ray detected. All 
but one of the PL classified AGN fall inside the 
Stern 'wedge' (Figure |4j, as expected from the 
non-independence of both indicato rs (as can be 



seen in the detailed comparison by iDonlev et al 



2008) 



The high AGN fraction and its increase with 
redshift, might lead one to think that the MIR 
criteria are overestimating the number of AGN at 
high redshifts, even though a tentative correction 
was applied (see Section f3.3.3|) . We have investi- 
gated the AGN fraction evolution from 1 < z < 2 
to 2 < z < 3 based, independently, on the X- 
ray and MIR indicators. In both wavebands, the 
AGN fraction increases significantly from low to 
high redshifts, rising from 9% to 16% when the 
X-ray is considered and from 15% to 29% when 
the MIR is considered. Although it seems possi- 
ble that SF galaxies may still be affecting the MIR 
criteria at high redshift (see Section 13.3. 3|) , most 
of the AGN fraction increase is consistent between 
the two wavelength indicators. This may partly be 
an effect of Malmquist bias, with lower luminos- 
ity systems, more likely to be dominated by SF, 
being progressively lost at higher redshift. In any 
case, this increase is consistent with t he known his 



tory of AGN activity in the universe ([Shaver et al 



19961 iHopkins et al. Il2007h 



The AGN fraction is not a simple function 
of color, as shown in Figure [H All three sub- 
populations (EROs, IEROs, and DRGs) have sim- 
ilar AGN fractions around the color threshold (~ 
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30%), but show different behavior with increasing 
colors. The more extreme colors among the ERGs 
do not necessarily correspond to higher fraction of 
AGN identifications, and in fact, the opposite ap- 
pears to be true for EROs. The AGN fraction is 
high only for the more "moderate" i — K s colors, 
decreasing for i — K s > 3.7. The difference be- 
tween the ERO and DRG trends result from each 
criterion itself. Redder i — K s colors will always 
select a population mostly at low- z (1 < z < 2), 
where the AGN fraction is shown to be smaller. 
On the other hand, redder J — K s constraints im- 
ply a higher fraction of high-z (2 < z < 3) sources, 
where the AGN fraction is higher. For example, 
essentially no low- z source has J — K s > 1.8 (Fig- 
ure!!)]). The AGN fraction versus color trend of the 
IEROs lies between that of the EROs and DRGs, 
which could be driven by the redshift distribution 
of the IERO population, which also lies between 
that of the EROs and DRGs. 

4.3. Radio Stacking 

Another important aspect necessary in under- 
standing the properties of ERGs is their SFR and 
the contribution of these populations to the overall 
SFR density of the universe. Dust obscuration is 
a serious source of uncertainty in estimating SFRs 
in ERGs from rest-frame ultraviolet luminosities. 
A SF diagnostic not affected by dust obscuration 
is radio emission. However, these galaxies are dis- 
tant enough that even the deepest radio surveys 
are sensitive to only the brightest star-forming 
population (detection limits corresponding to sev- 
eral hundred M Q yr _1 for z > 1). Instead, stack- 
ing methods can be used to evaluate the statis- 
tical star-forming properties of ERGs. Stacking, 
as used here, is simply an "image stacking" pro- 
cedure, where image sections centered at each de- 
sired source position (stamps) are combined. The 
aim is to reach much lower noise levels, possibly 
providing a statistical detection of samples whose 
elements are individually undetected in the origi- 
nal image. 

For the radio stacking analysis we have used the 
1.4 GHz ATCA observations of this field, reaching 
a uniform 1 4-17 nJy rms throughout the GOODS- 
S field (see lAfonso et al.ll2006t iNorris et al]|2006l 
for more details). Our adopted stacking method- 
ology can be summarized in the following steps. 

First, using the radio image of the field, stamps 



of 60 by 60 pixels (equivalent to 120" by 120") 
were deemed appropriate, allowing a good sam- 
pling of the vicinity of each source, necessary to 
identify strong neighboring sources that can bias 
the stacking. 

Every stamp containing a radio source within 
an 18" radius from the central (ERG) position, 
was rejected. As the stacking will be used to es- 
timate the average flux of the unidentified ERGs 
in the radio image, the inclusion of radio detec- 
tions would likely bias the final result. In this 
context the term "detection" does not only apply 
to the robust detections (roughly at > 4.5a level), 
but also to the "possible" detections (all remaining 
candidate radio sources at > 3cr). Stamps having 
radio sources near the ERG position must also be 
excluded, as the wings of the radio detection can 
extend to the central part of the stamp. 

The remaining stamps for each sample of ERGs 
can then be stacked. Previous work often uses 
median stacking (e.g., 



White et al 



20071 ). in an 

attempt to be robust to radio detections and 
high/low pixels. The penalty for this is the loss 
of sensitivity. Having removed all detections and 
possible detections from the list of stamps, we use 
a weighted average (weight=rms -2 ) stacking pro- 
cedure. At each pixel position we further imple- 
ment a rejection for outliers, rejecting high (low) 
pixels above (below) the 3<r (—3c) value for that 
pixel position. The number of rejected pixels in the 
central region is always very low (< 2) confirming 
that previous rejection steps work efficiently. 

The final flux and the noise level are measured 
in the resulting stacked image. To evaluate the 
reliability of detections in the stacked images we 
performed Monte Carlo (MC) simulations. Ran- 
dom positions in the radio image were selected and 
stacked, following the procedure described above. 
Each of these positions were required to be far- 
ther than 6"from the known K s sources, as we 
are interested in evaluating systematics of the ra- 
dio image alone. Appropriate numbers of stacked 
stamps were used, to compare to the actual ERG 
stacks. The procedure was repeated 10000 times 
for a given number of stamps. For further refer- 
ence, a stacked sample will be considered to have 
produced a reliable detection only if no MC sim- 
ulation (among 10,000) has resulted in a higher 
S/N value. 
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4.4. Star formation activity in ERGs 

Following the procedure outlined above, we 
have performed a radio stacking analysis for differ- 
ent sub-groups within the ERG population. The 
radio data were stacked for each of the popula- 
tions of EROs, IEROs, DRGs, pEROs, pDRGs, 
and cERGs. Within these samples, stacking of 
the radio images was also performed separately for 
AGN and non-AGN sub-populations. Since red- 
shift estimates exist for the vast majority of the 
ERGs, stacking is performed separately for both 
low and high redshifts (1 ^ z < 2 and 2 ^ z ^3, 
respectively). Besides minimizing biases in the 
stacking signal, due to different populations and 
different (radio) luminosities being sampled at dif- 
ferent redshifts, this also allows us to search for a 
hint of any evolutionary trend. Given the incom- 
pleteness of the sample at the highest redshifts 
no attempt was made to perform a specific radio 
stacking analysis for z > 3 ERGs. Table [2] lists 
the number of sources considered in each of the 
sub-populations and those in each of the stacking 
steps referred in the previous section. 

While the stacking procedure enables the aver- 
age flux to be estimated from the radio-undetected 
sample (< 3cr signal), the entire population should 
be considered when measuring the ERG contribu- 
tion to the global SFR density of the universe. The 
approach adopted here was to consider all radio- 
undetected ERGs as having a radio flux given by 
the average signal from the stacking analysis, and 
all radio-detected ERGsJjl to contribute with their 
measured flux density. The conversion from radio 
flux to radio luminosity is performed by using the 
assumed redshift (spectroscopic or photometric) 
and a radio spectral index of a — 0.8 (S v oc v~ a ): 

Lxaguz = 4tt4 Si.4ghz 10~ 33 (1 + z)"- 1 WHz-\ 

where is the luminosity distance (cm) and 
»5i.4GHz is the 1.4 GHz flux density (mJy). The 
correspon ding S FR is obtained using the calibra- 



tor individual galaxies using the 1/V ma x method 



(Schmidt 1968). The associated volume for each 



tion from [Belli ([20031 ): 



SFR (M Q yr" 1 ) = 



5.52 x 10"" 

-1-22 



5.52x10' 



I 0.1+0.9(^) 

where L c = 6.4 x 10 21 W Hz" 1 
Hz -1 . The contribution to 



il.4GHz 
trail. 4GHz 



L> L c 
L<L r 



= 10 21 - 81 W 
was estimated 



galaxy is estimated by using a i^-correction de- 
rived from the galaxy's own SED (as given by the 
observed multi- wavelength photometry). Again, 
radio-detected ERGs, contributed with their esti- 
mated intrinsic luminosity and SFR, derived with 
the assigned redshift estimate and its detected 
flux. In Figure [101 the SFR distribution of these 
sources is presented. Those classified as SF (14 
in total), range from ^50 to ^1600 M Q yr -1 
(in reasonable agreemen t with those presented in 
Georgakakis ct al. 2006). On the other hand, the 



luminosity and SFR estimates of radio- undetected 
ERGs were based on the resulting stacking signal 
of the sample and, likewise, the individual ERG 
redshift value. 

The results are given in Table [3] For each 
ERG sub-population we list: (1) the ERG sub- 
population; (2) the total number of sources in the 
sample; (3) the final number of stamps included 
in the stacking; (4) the rms of the final stacked 
image; (5) the measured flux in the central region 
of the stacked image; (6) the respective S/N; (7) 
number of Monte-Carlo simulations (out of 10000) 
that resulted in higher S/N values, a measure of 
the reliability of the ERG detection: conserva- 
tively, whenever Nmc > the stacking signal is 
considered spurious; (8) the average redshift for 
the sub-population; (9) the average radio luminos- 
ity both for the radio non-detected sources - tak- 
ing into account the stacking signal only - and, in 
parenthesis, that for the entire sub-population (in- 
cluding radio detected sources); (10) the average 
SFR, for non-AGN samples, corresponding to the 
radio luminosities in column (9); (11) the result- 
ing radio luminosity density (£1.4 GHz); (12) the 
corresponding p*. For columns (9) to (12), the 
upper limits, corresponding to non-detections of 
the stacked signal (Nmc > 0), are estimated using 
the maximum S/N found on the MC simulations. 
No stacking is attempted for populations with less 
than 10 stamps. 

The analysis suggests that the bulk of the ERO 
populations have modest SF activity. At 1 < z < 
2, where most EROs are found, the average SFR 



4 For this purpose, radio-detections refer to signals above 3a 
in the radio map; see Section 14.31 



5 In this section, while discussing the SF properties of ERGs, 
we will naturally be referring to the non-AGN sources, 
dropping the prefix 'n' for simplicity. 
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is below a few M© yr" 1 . Only at 2 < z < 3, 
EROs - practically all (93%) being classified as 
DRGs - reveal intense average SFRs, ~ 200 Mq 
yr -1 , entering the Ultra Luminous Infrared Galax- 
ies (ULIRG) regime. This suggests that at \ow-z 
the passive/evolved systems represent a significant 
fraction of the ERO population (56%, see pEROs 
discussion ahead) , as opposed to the high-z regime 
where the dusty systems dominate. DRGs show 
the same trend although not as pronounced: they 
increase from ~ 70 Mq yr -1 at 1 < z < 2 to ~ 140 
Mq yr -1 at 2 < z < 3. Interestingly, at low red- 
shifts (1 < z < 2), 97% of DRGs (69 out of 71) are 
also classified as EROs; conversely only 20% of the 
low-z EROs (69 out of 344) are also classified as 
DRGs. The estimated average SFR of ~ 70 M 
yr" 1 for this overlapping population supports pre- 
vious cl aims of a dusty starburst nature for the se 
sources (jSmail et al.ll2002l : IPapovich et al.ll2006l ). 

These values for the DRG population arc 
comparable to those found in the literature 
dRubin et al.l l2004UForster Schreiber et all 12004 



Knu dsen et al .l20a4lR,eddv et al.l2 005) . IPapovich et al 
(|2006l ) studied 153 DRGs selected also in the 
GOODS-S to a limiting magnitude of K Si tot < 
23 AB. They find an average SFR for the DRG 
population at 1 < z < 3 of 200 - 400 Mq yr" 1 , 
which is somewhat higher than our result. This 
difference is attributable to the different meth- 
ods of DRG selection. Firs t, the DRG selec- 
tion in IPapovich et al. (l2006h is slightly redder 
(0.06 mag) and restricted to brighter if-band 
magnitudes (0.6 mag difference). The first alone 
would be responsible for an increase of ~10% in 
our estimated SFR, while the latter effectively 
restricts the sample to stronger starbursts, in- 
creasing the average SFR by ~30%: from 100 
M Q yr" 1 at K s ^ TO t < 23. 8 to 130 M m yr" 1 



at K s ^ T ot < 23.0(see also iReddv et ail (20051 ). 
Second, we employ a more conservative AGN re- 
jection, by using the full 2Ms Chandra data and 
rejecting essentially all sources cl assified as AGN 
from MIR color-color diagrams. IPapovich et al 



(|2006l ) look at the individual MIR IRAC bands to 
re-classify as starforming a considerabl e numb er of 
DRGs flagged as AGN bv lStern et all (|2005h . Fi- 
nally, the uncertainty in the different methods and 
associated calibrations (UV+IR luminosities ver- 
sus Radio luminosities) used to estimate the SFR 
account for the remaining observed difference. 



The low average SFR for EROs at 1 < z < 2 
is due to the numerous pEROs (193, 56% of the 
1 < z < 2 EROs) : the stacking analysis of pEROs 
found in this redshift range fails to produce any 
signal. This population likely corresponds to the 
passively evolving component of EROs. On the 
other hand, pDRGs at 2 < z < 3 must also be 
characterized by relatively low SFRs: although the 
stacking analysis is unable to give such indication 
(only constraining the average SFR to < 300 Mq 
yr -1 ), pDRGs are the sources responsible for the 
observed difference of the average SFR of EROs 
and that of DRGs in this redshift range (all but 
4 of these non-AGN EROs are also classified as 
DRGs). The pDRGs effectively dilute the stack- 
ing signal from 160 (200 if radio detections are 
included) M Q yr" 1 to 110(140) M Q yr" 1 . The 
latter should thus be taken as a better upper SFR 
limit than that obtained from the actual pDRG 
stack. 

The SFR density behavior for ERGs roughly 
follows the general trend for SF galaxies, increas- 
ing from l<z<2to2<z<3 (Figure fTTj) . Over- 
all, the ERG contribution to the total p* jumps 
from ~ 10% in the low redshift bin to around 
20% at 2 < z < 3, where EROs are the highest 
contributors (~ 0.043 Mq yr" 1 Mpc" 3 ). The ex- 
ception is the IERO population which apparently, 
despite being the biggest contributor among the 
three populations at low-z, shows no evolution be- 
tween the two redshift intervals. 

We finally highlight the fact that these esti- 
mates may be missing potentially significant SF 
from the ERGs classified as AGN. The imple- 
mented AGN identification criteria were applied 
in a very conservative way, possibly even removing 
some non-AGN galaxies from the non-AGN sam- 
ples. Also, the SFRs associated with the AGN 
populations are obviously being discarded in this 
as in similar investigations. This separation is 
currently necessary given the lack of knowledge 
on how to remove the AGN contribution to the 
(radio) emission. The lack of powerful AGN in 
the ERG populations suggest omitting any at- 
tempted AGN exclusion, without significant detri- 
ment to the average SFR estimates: for example, 
the entire ERO population in the 2 < z < 3 red- 
shift range displays an average radio luminosity of 
10 23 51 WHz" 1 , only 10% higher than that for the 
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non-AGN EROs in the same redshift ranged. How- 
ever, when no AGN discrimination is attempted, 
the numerous AGN sources would result in a much 
more significant increase for the SFR density es- 
timate: in the previous example, the SFR density 
for the ERO population would rise by > 100%, 
from 0.04 to 0.1 M yr -1 Mpc -3 . The correct sce- 
nario will obviously be between the two extremes. 
A precise handling of such uncertainty can only be 
fully addressed once a significantly better under- 
standing of the AGN-SF connection is achieved, 
and mor e discriminating dia gnostics are developed 
(but see lDunne et al. 1 120091) . 



5. Conclusions 

We have presented a multi- wavelength analysis 
of the properties of the ERG population in the 
GOODS-South field. EROs, IEROs, DRGs - and 
various combinations between these groups - are 
considered, their AGN content identified and their 
contribution to the global p» estimated, leading to 
the following conclusions. 

• The different criteria for the selection of red 
galaxies select, as previously known, sources 
at different redshift ranges: while the bulk of 
EROs and IEROs can be found at 1 < z < 2, 
DRGs are mostly found at 2 < z < 3. Differ- 
ent combinations of the three criteria result 
in samples with distinct redshift properties: 
while cERGs are observed in a wide redshift 
range, 1 < z < 3, and have no low-z (z < 1) 
interlopers, pEROs and pDRGs appear in 
distinct redshift intervals, at 1 < z < 2 and 
2 < z < 4, respectively. The "pure" criteria 
appear, thus, to be suitable and simple tech- 
niques to select high-z sources in well con- 
strained redshift intervals. See Section Hj] 



the ERG population does not include a large 
number of powerful AGN, as indicated by 
the X-rays and radio observations. Almost 
one-third of the ERG sample hosts poten- 
tial AGN activity, with the fraction of AGN 
increasing from EROs to IEROs to DRGs 
(respectively, 29%, 35%, and 38%). Among 
ERGs, and according to the X-ray proper- 
ties, type-2 sources dominate (a 2:1 ratio). 



An X-ray estimate of the type-2 to type- 
1 AGN ratio among the ERG population 
is, however, indeterminate, requiring obser- 
vations extending to lower X-ray energies 
(higher wavelengths). See Section 13.21 and 
Section l4~2l 



• EROs at z < 2 are often pEROs (56%), 
which are likely passively evolved systems 
without strong SFR activity, on average be- 
low ~ 20 Mq yr -1 . On the other hand, es- 
sentially all EROs at 2 < z < 3 are classified 
as DRGs and have intense SF activity, simi- 
lar to ULIRGs (on average ~ 200 M Q yr^ 1 ). 
DRGs also show this increase in average SFR 
from low (1 < z < 2) to high (2 < z < 3) 
redshifts, although not as pronounced (from 
- 70 to - 140 M© yr" 1 ). See SectionSU 

• The overlapping population between EROs 
and DRGs displays an intense average SFR 
both at 1 < z < 2 (~ 70 M yr" 1 , sup- 
porting previous claims o f a dusty starburst 



nature for these sou rces, ISmail et al. 2002; 



Papovich et al J 120061 ) and 2 < z < 3 (~ 200 



M Q yr" 1 ). See Section HH 

• The contribution of ERGs to the SFR den- 
sity increases with redshift: from ~ 10% at 
1 < z < 2 to ~ 20% at 2 < z < 3, an 
increase which is attributable to the ERO 
and DRG populations. IEROs, on the other 
hand, despite showing the highest contribu- 
tion to the global SF history among the three 
ERG population at low-z, appear to show no 
evolution at high redshifts. See Section FOl 

• SFR densities from ERG populations were 
estimated for SF-dominated sources, after 
a thorough AGN multi-wavelength identifi- 
cation. We find that, although the AGN 
ERGs would only slightly increase the av- 
erage radio luminosities shown by the non- 
AGN samples, inclusion of such sources in 
the SFR density estimates could lead to sig- 
nificant (and, at this point, undetermined) 
biases. See Section [4~4l 



6 Note, however, that the increase for IEROs at 1 < z < 2 
reaches > 30%. 



The use of a [5.8]-[8.0] versus [8.0]- [24] color 
diagnostic allows for a tentative separation 
between AGN and SF galaxies at z > 2.5, 
where mid-IR (3 to 8 /im) diagnostics be- 
come degenerate. In particular, the use of 
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this diagnostic enables the identification of 
a z ~ 2.5 evolved system candidate with 
MIR colors typical of Sc galaxies. See Sec- 
tion [3331 
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DRGs 




Fig. 1. — Adopted ERG sub-sample nomenclature 
is represented here through a Venn diagram. The 
overlap between the three ERG classes - EROs, 
IEROs, and DRGs - is significant (the common 
ERG population, labeled as cERGs). The outer 
non-overlapping regions represent the pure popu- 
lations. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 2. — K s — [3.6] color-magnitude plot for 
sources in the FIREWORKS catalog. Points and 
open histograms (scaled down by a factor of 4) rep- 
resent the general K s population included in the 
catalog, while filled circles and shaded histograms 
represent the selected IERO sample. The diagonal 
line corresponds to a i^-band value of 23.8 mag, 
the adopted magnitude limit of our sample. The 
dashed vertical line represents the practical [3.6]- 
band IERO magnitude limit (as imposed by the 
Z850 3<7 magnitude limit and the IERO definition). 
The current AVdctcctcd IERO sample would dif- 
fer significantly from the general IERO population 
if a large number of sources exists above the diag- 
onal and to the left of the vertical line (shaded 
region). The K s — [3.6] colors present in the well 
sampled region of the diagram (below the diago- 
nal line and to the left of the vertical one) argue 
against this scenario, implying that the current 
sample of IEROs is representative of the overall 
IERO population. 
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Fig. 3. — X-ray HR evolution with redshift for ob- 
scured (Nji = 10 23 cm -2 , gray shaded region) and 
unobscured (Njj — 10 2CK2 cm~ 2 , light gray shaded 
region) X-ray PL emission models (r = 1.8 ±0.5), 
calculated using PIMMS (ver. 3.9k). Filled cir- 
cles show the distribution of the X-ray detected 
AGN ERGs with a robust HR estimate. Upper 
limits (no hard-band detection) appear as empty 
triangles while filled triangles denote lower limits 
(no soft-band detection). The dashed horizontal 
line highlights the HR const raint (HR= —0-2) fo r 
type discrimination used bv lSzokolv etaD (|2004h . 
It is clear that for high-redshift sources (z > 2) 
the simple HR criterion becomes degenerate as an 
obscuration measure. 
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. ERGs 
# PL AGN 

A XR AGN 



[5.8]-[8.0] (AB) 

Fig. 4. — Distribution of ERGs on the colo r -color 
diagnostic plot proposed by IStern et ah ( 2005 ) 



with the AGN region delimited by the dashed 
lin e. The ERGs c lassified as AG N by the PL 



I Do nley et al.l 120071 ) and X-rays (lSzokolv et al 



2004) criteria are highlighted (red circles and green 
triangles, respectively). 



[5.8]-[8.0] 

Fig. 5.— Mid-infrared [5.8]- [8.0] versus [8.0]- [24] 
color-color plot for z > 2.5 ERGs in the cur- 
rent sample. Filled symbols represent an AGN 
classification from a MIR indicator (IRAC PL or 
IRAC color-color) and open symbols otherwise. 
The tracks represent the expected colors of tem- 
plate SEDs where the IR is dominated by star- 
formation (dotted and dashed tracks, where the 
latter represent more intense SF activity) or AGN 
activity (continuous tracks), redshifted between 
z = 2.5 and z = 4, with crosses at z = 2.5 
and z = 3. The templates displayed are (from 
bottom to top): 2 Spirals (Sc and Sd), 3 star- 
bursts (M82, NGC 6240, and Arp220), and 6 
A GN (IRAS 22491-18 08, IRAS 20551-4250, QSO- 
2, lAfonso et all (|200lh ERO, Mrk 231, and IRAS 
19254-7245 South). 
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Fig. 6. — Redshift distributions of different ERG 
sub-populations: (a) EROs, (b) IEROs, and (c) 
DRGs. The hatched histograms correspond to 
ERGs identified as AGN. Note the different y-axis 
scales for the individual panels. 
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Fig. 8. — AGN fraction as a function of color for 
EROs, IEROs, and DRGs. The x-axis represents 
the difference in color to the color threshold ad- 
equate for each population: 2775 — K s = 2.48 for 
EROs, x 850 - [3.6] = 3.25 for IEROs and J-K s = 
1.3 for DRGs. IEROs appear to have an inter- 
mediate behavior relative to EROs and DRGs, as 
one goes to more extreme colors. The histograms 
on the upper part of the picture show the color 
distribution for each ERG sub-population. 



Fig. 7. — Redshift distributions for "pure" and 
"common" ERG sub-populations: (a) pEROs, (b) 
pDRGs and (c) cERGs. The hatched histograms 
correspond to ERGs identified as AGN. Note the 
different y-axis scales for the individual panels. 
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.AGN DRGs 

o non-AGN DRGs 




Fig. 9. — Variation of J — K s color with redshift. 
The DRG criterion color cut is shown as a dashed 
line. A 0.5 mag redder J — K s cut (dotted line) 
selects almost no z < 2 DRGs. 




Z 6 - 




Fig. 11. — Contribution of ERG populations to 
the total p» for 1 < z < 2 and 2 < z < 3. EROs 
are denoted by green diamonds, IEROs by cyan 
triangl es, and DRGs by blue circle s. The compila- 
tion of iHopkins fc Beacom I (|2006f) is displayed for 
reference (gray crosses and shaded region, corre- 
sponding to the p* la and 3er confidence regions). 



sfr i.4ch« [ M o yr" 1 ] 



Fig. 10. — SFR distribution of the radio detected 
ERGs (considering any signal in the radio map 
with > 3a). Dashed histograms show the overall 
distribution, while continuous histograms refer to 
the sources considered as SF systems. 
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Table 1 

ERG Number Statistics: Sample Overlap, Counterparts, and Classification. 



POP ATtot 3 ero iero drg JV 3pcc X-Ray IRAC d PL C MIR cc° MIPS Radio N 

XR b Al c A2 C Ql c Q2 C nHR c 24/1111 1.4 GHz 



ERG 


731 (702) 


645 


294 


350 


155 


88 


18 


36 


1 


5 


26 


675 


14 (1) 


164 (30) 


402 


17 


214 (29%) 


ERO 


645 (626) 


645 


290 


267 


145 


84 


17 


35 


1 


4 


25 


612 


6 (1) 


135 (27) 


359 


15 


184 (29%) 


IERO 


294 (289) 


290 


294 


198 


29 


47 


9 


19 


1 


3 


14 


290 


4 (1) 


83 (21) 


191 


10 


103 (35%) 


DRG 


350 (338) 


267 


198 


350 


39 


51 


11 


17 





5 


16 


322 


12 (0) 


117 (25) 


222 


11 


135 (39%) 


cERG 


197 (197) 


197 


197 


197 


17 


36 


8 


14 





3 


10 


197 


2 (0) 


69 (18) 


141 


7 


81 (41%) 


pERO 


272 (272) 


272 








103 


24 


6 


13 








5 


257 


(0) 


34 (2) 


128 


3 


56 (21%) 


pDRG 


70 (70) 








70 


9 


4 


1 


1 





1 


1 


51 


6 (0) 


24 (3) 


35 


2 


25 (36%) 



a Total number of sources in each (sub)samplc and, in parenthesis, those which have good photometry in all bands involved in the ERG criteria: ^775, zq^q. J. K s , 
and 3.6p.m. 

b Total number of X-ray identifications. 

c Number of sources classified as typc-1 or type-2 A GN (Al or A2, respectively), type-1 or typc-2 QSO (Ql or Q2, respectively), and AGN with undetermined type 
(no HR determination, column nHR) according to the Szokoly ct al. (2004) criterion. 

d Number of sources with detection in all four IRAC channels. 

c Sources classified as AGN by PL (Donley ct al. 2007) or MIR color-color (Stern ct al. 2005) criteria. In parenthesis, those also classified as AGN by the X-ray 
criterion. 

*Total number of sources classified as AGN, considering all AGN identification criteria, along with the equivalent fraction in the total (sub)population. 



Table 2 

Robust Radio stacking of ERG populations. 



POP Wtot ^V<i8'' a ^3<, b JV fi „ c 



EROs 










zl2 


443 


56 


12 


370 


zl2; nAGN 


344 


37 


7 


296 


z23 


144 


16 


3 


124 


z23; nAGN 


84 


6 


2 


75 


IEROs 










zl2 


161 


20 


6 


132 


zl2; nAGN 


124 


11 


5 


106 


z23 


111 


12 


2 


96 


z23; nAGN 


59 


3 


1 


54 


DRGs 










zl2 


102 


9 


3 


88 


zl2; nAGN 


71 


4 


2 


64 


z23 


176 


20 


3 


152 


z23; nAGN 


111 


7 


2 


101 


cERGs 










zl2 


71 


7 


2 


60 


zl2; nAGN 


49 


3 


2 


43 


z23 


104 


12 


2 


89 


z23; nAGN 


59 


3 


1 


54 


pEROs 










zl2 


245 


34 


5 


204 


zl2; nAGN 


193 


25 


2 


164 


z23 


5 


1 





4 


z23; nAGN 


4 


1 





3 


pDRGs 










zl2 


1 











zl2; nAGN 


1 











z23 


34 


4 





30 


z23; nAGN 


25 


1 





24 



Note. — The zl2 and z23 abbreviations stand for 
1 < z < 2 and 2 < z < 3, respectively. 

a Number of stamps with a radio detection within 
18" of the ERG position, consequently rejected from 
the final stacking. 

b Number of stamps with a possible radio detection 
at the ERG position (signal between 3<r and ~ 4.5<r), 
also removed from the final stacking. 

c Final number of stamps included in the stacking. 
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Table 3 

Properties of Extremely Red Galaxy Populations from Radio Stacking Analysis 



to 



POP 


JVtot 


Affi„ a 


rms 


Sl.4GHz 


S/N 


N MC h 


z 


log (Ll.4GHz) C 


SFR C 


log (£l.4GHz) C 










(z^Jy) 


!>Jy) 








(\VHz~ 1 ) 


(Moyr- 1 ) 


(WHz -1 Mpc~ 3 ) 


(Afgyr^Mpc" 3 ) 


EROs 
























zl2 


443 


370 


0.927 


2.051 


2.212 





1.45 


22.37(22.83) 




19.41(19.90) 




zl2; nAGN 


344 


206 


1.037 


1.080 


1.041 


5 


1.45 


<22. 25(22. 62) 


<10(23) 


<19. 19(19. 56) 


<8.48e-03(2.01e-02) 


z23 


144 


124 


1.708 


8.641 


5.059 





2.44 


23.51(23.62) 




20.05(20.17) 




z23; nAGN 


84 


75 


2.006 


8.274 


4.125 





2.39 


23.47(23.55) 


162(196) 


19.81(19.90) 


3.60e-02(4.34o-02) 


IEROs 
























zl2 


161 


132 


1.627 


6.700 


4.118 





1.66 


23.01(23.25) 




19.63(19.96) 




zl2; nAGN 


124 


100 


1.792 


5.071 


2.830 





1.67 


22.89(23.08) 


43(66) 


19.42(19.62) 


1.45e-02(2.32o-02) 


z23 


111 


90 


1.928 


8.176 


4.241 





2.47 


23.50(23.58) 




19.94(20.03) 




z23; nAGN 


59 


54 


2.369 


7.103 


2.998 


1 


2.40 


<23. 46(23. 49) 


<161(170) 


<19. 68(19. 71) 


<2.67e-02(2.82o-02) 


DRGs 
























zl2 


102 


88 


2.001 


6.823 


3.410 





1.59 


22.99(23.12) 




19.39(19.56) 




zl2; nAGN 


71 


64 


2.283 


6.711 


2.940 





1.61 


22.99(23.09) 


54(68) 


19.27(19.43) 


1.04e-02(1.48e-02) 


z23 


17G 


lo2 


1.596 


6.144 


3.850 





2.51 


23.39(23.52) 




19.95(20.10) 




z23; nAGN 


111 


101 


1.744 


5.224 


2.995 





2.47 


23.30(23.40) 


110(137) 


19.69(19.79) 


2.73e-02(3.42e-02) 


cERGs 
























zl2 


71 


GO 


2.542 


8.901 


3.502 





1.67 


23.15(23.26) 




19.49(19.64) 




zl2; nAGN 


49 


43 


2.982 


7.275 


2.440 





1.71 


<23. 16(23. 26) 


<80(100) 


<19. 39(19. 55) 


<1.36e-02(1.96o-02) 


z23 


104 


89 


1.961 


7.270 


3.707 





2.48 


23.45(23.55) 




19.86(19.97) 




z23; nAGN 


59 


54 


2.369 


7.103 


2.998 


1 


2.40 


<23. 46(23. 49) 


<161(170) 


<19. 68(19. 71) 


<2.67e-02(2.82e-02) 


pEROs 
























zl2 


245 


204 


1.289 


-0.449 


-0.348 


1347 


1.33 


<22. 36(22. 66) 




<19. 17(19.47) 




zl2; nAGN 


193 
6 


164 


1.413 


-1.036 


-0.733 


3764 


1.32 


<22. 52(22. 68) 


<18(26) 


<19. 24(19. 39) 


<9.54e-03(1.36e-02) 


z23 
z23; nAGN 


4 






















pDRGs 
























zl2 


1 






















zl2; nAGN 


1 






















z23 


34 


30 


3.713 


1.049 


0.283 


3396 


2.67 


<23. 71(23. 75) 




<19. 67(19. 72) 




z23; nAGN 


25 


24 


4.037 


-1.532 


-0.379 


7001 


2.65 


<23. 80(23. 80) 


<345(345) 


<19. 70(19. 70) 


<2.77c-02(2.77e-02) 



Note. — The zl2 and z23 abbreviations stand for 1 < z < 2 and 2 < z < 3, respectively. The upper limits for Luminosity and SFR estimates, whenever 
Nmc > 0, arc calculated considering the maximum S/N obtained in the respective set of MC simulations. 

a Final number of stamps included in the stacking, after the various rejection steps described in Scction |4.3l 

b Number of MC simulations (out of 10000) that resulted in higher S/N values. 

L In parenthesis, the estimated value also takes into account radio detections (> 3(t) excluded from the stacking procedure (see Section I4.3|) . 



